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From study of the Kosterlitz-Thouless-Berezinskii (KTB) transition in the superfluid density, 
n B (T), of ultrathin c-axis oriented YBa2Cu307_a (YBCO) films, we find that interlayer coupling is 
unexpectedly strong. The KTB transition occurs at a high temperature, as if the films were isotropic 
rather than quasi-two-dimensional. This result agrees with a comparison of the superfluid density 
of YBCO with Bi2Sr2CaCu2 08 and with numerical simulations of Josephson junction arrays, and 
challenges the thermal phase fluctuation interpretation of critical behavior near T c in YBCO. 

PACS numbers: 74.25.Fy, 74.40. +k, 74.76. Bz, 74.72.Bk 
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Studies of thermal fluctuations in the phase of the superconducting order parameter have a long history due to 
their importance in understanding phase transitions in general. Because of the high superfluid densities, n s , and long 
coherence lengths, £, of conventional superconductors, fluctuations in three-dimensional (3D) samples are significant 
only within a hopelessly narrow temperature range near T c p| , and most experimental studies involve thin amorphous 
films in which reduced dimensionality and superfluid density enhance fluctuations to an observable level. Two points 
that are relevant to work on fluctuations in high-T c cuprate superconductors have been established. First, quantum 
mechanical effects suppress thermal phase fluctuations (TPF's) when T drops below a surprisingly high temperature 
[2-|4|. Second, at high temperatures where fluctuations are classical, simulations of Josephson junction (J J) arrays 
[5| are in reasonable quantitative agreement with measurements of the superfluid density in two-dimensional (2D) 

vq . films J2| . In this paper we concern ourselves with fluctuations in hole-doped cuprate superconductors in zero applied 

f^i ' magnetic field. 

Thermal phase fluctuations were immediately expected to be significant in cuprate superconductors |L| , which have 
low superfluid densities and weak interlayer coupling. Exactly how important TPF's are is a matter of considerable 
debate. In the absence of directly comparable experimental results on quasi-2D conventional superconductors, we rely 
on simulations of strongly anisotropic JJ arrays p] to provide a framework for examining experimental results. In 
YBa2Cu307_a (YBCO), strong fluctuations are evidenced by very rapid changes in the a&-plane penetration depth 
X(T) J7H9], specific heat ]l0| , and coefficient of thermal expansion |1J] over a 5 to 10 K interval near T c . Perhaps 
the most striking evidence that these phenomena are due to TPF's is power-law behavior with the 3D-XY exponent 
2/3 in the superfluid density, n s (T) ex A _2 (T), of YBCO crystals |7]-||. In the TPF interpretation, fluctuations are 
strong because of quasi-two-dimensionality, and interlayer coupling - although weak - begets a crossover from 2D to 
the observed 3D fluctuations |qjl2|. 

Why not just accept these results at face value and move on? The reason is that T c marks a second order phase 
transition from an unusual normal state to a d-wave superconductor. It would not be surprising if normal-state 
peculiarities persisted into the superconducting state and resulted in anomalous T-dependence of n s . For example, if 
there is a normal-state order parameter (e.g., see Ref. 13), it may compete with the superconducting order parameter 



In the present paper, we argue from measurements of n s in the literature and new measurements presented below, 
that TPF effects are actually small, presumably because interlayer coupling in cuprates is much stronger than indicated 
by transport measurements. 

The pivotal notions of the TPF model are that fluctuations are strong because interlayer coupling is very weak, and 
that interlayer coupling serves only to change the T-dependence of n s (T) from a discontinuous Kosterlitz-Thouless- 
Berezinskii (KTB) |]15| drop at T^tb t° a continuous decrease to zero, with rapidly increasing downward curvature 
p,[L2| . One key quantity in comparing n s (T) with the TPF model is the hypothetical transition temperature, T^ TB , 
at which the KTB transition would occur if layers could be completely decoupled without affecting the pairing 



interaction. T^ TB can be deduced from the measured a6-plane superfluid density. A quantitative measure of the 
strength of fluctuations, and therefore of interlayer coupling, is the difference between T^ TB and the temperature, 
T c , at which n s (T) vanishes: the weaker the coupling, the smaller the difference. From numerical results in Ref. 6, 
T c /Tf CTB — l increases roughly as (Jj_/J||) 1/ ' 3 , where J± (Jm) is the interlayer (intralayer) coupling strength. The other 
key quantity is n s (T^ TB ) oc A~ 2 (Tj*- TB ). In the TPF model, 3D-XY critical behavior develops only as n s approaches 
zero, i.e., for n s (T) <C n s (T^ TB ). 

Let us now take a closer look at power-law behavior in n s {T) observed for some YBCO samples, setting aside 
concerns that power-law behavior is not seen in other YBCO samples |l6| , |l7| ] nor in clean Bi2Sr2CaCu20s (BSCCO) 
crystals ]18[ ]. A recent measurement of the a-axis magnetic penetration depth in a clean, optimally doped YBCO 
crystal [Bj finds that at temperatures from ~ 85% to ~ 99.9% of the transition at T c rs 93.78 K, X~ 2 (T) is fitted to 
within a percent or so by A~ 2 (T)/A~ 2 (0) rs 1.26(1 — T/T c ) 2 / 3 . The critical region encompasses a remarkably large 
range of temperature and superfluid density: 0.85 < T/T c < 0.999 and 0.01 < A~ 2 (T) / A~ 2 (0) < 0.35. The argument 
that an excellent fit is unlikely to be accidental is taken as support for the TPF model, but in reality it is not. The 
reason is that in the TPF model the crossover from 3D to 2D fluctuations occurs well inside the observed "critical 
region". To see this, let us estimate T^ TB and A~ 2 (Tj^ TB ) from the data. With an effective layer thickness equal to 
the center-to-ccntcr spacing between Cu02 bilayers (d — 1.17 nm), A a (0) = 160 nm, and the KTB prediction |15[ 

A :'< T ™» = ss^p (1 » 

where AJ 1 = dA~ 2 , we find T* KTB w 0.96T C and K 2 ( t ktb) ~ ^a 2 (0)/7. The TPF model predicts 3D-XY behavior 
only for A~ 2 (T) < A~ 2 (T£ TB ), which for these data means T > 0.99T C . Power-law behavior from 0.85 to 0.99 T c 
would have to be "accidental". Taking this one step further, if power-law behavior is due to physics beyond the TPF 
model, then TPF's must be so weak as to be unapparent in the data. 

Another argument against the TPF model arises from comparison of YBCO with BSCCO, which has Cu02 bilayers 
that are nearly identical structurally to those of YBCO - i.e., the same interplanar spacing and nearly square Cu-0 
structure - except that Y 3+ is replaced by Ca 2+ p9| . As expected, the areal superfluid density, n s (T)d, at T = of 
the bilayers in clean BSCCO is essentially the same as in the a-axis direction (perpendicular to CuO chains) in YBCO. 
However, BSCCO has interlayer coupling several hundred times weaker than YBCO, as measured by transport pTj| , pT| , 
so we expect TPF's to be much stronger. Let us see. T c lies 4 K above T^ TB in YBCO crystals, as found above - a 
reasonable result when compared with simulations of JJ arrays 0. In BSCCO, the difference is 10 K, indicating that 
fluctuations in BSCCO are actually weaker than in YBCO. The latter number comes from clean BSCCO crystals 111] 
in which X~ 2 (T) decreases linearly with T at low T, T c is 91 K, the center-to-center spacing between Cu02 bilayers 
is 1.54 nm |y|, and A(0) = 210 nm. Evidently interlayer coupling is stronger in BSCCO than indicated by transport 
measurements of a6-plane vs. c-axis anisotropy. 

We could refine the above estimations by using two interlayer coupling constants, one within bilayers and another 
between bilayers, per Carlson et al. ||, and use an average of the a&-plane superfluid density rather than just the 
a-axis value in YBCO, but the conclusions do not change. 

Other experimental evidence supports the proposition that interlayer coupling is stronger than expected. Some 
corroboration is provided by a study of <ji(w,T) above T c in YBCO crystals ||, which finds a critical region of only 
0.6 K. More direct support is provided by the measurements presented below, which demonstrate that fluctuation 
effects in the 2D penetration depth, Aj/(T) oc ns{T)d, of ultrathin twinned YBCO films (d = 4, 8, and 10 unit 
cells) are quantitatively consistent with measurements on 2D s-wave films [g,^2[, assuming that the YBCO films are 
effectively isotropic - i.e., that their layers are strongly coupled. 

II. EXPERIMENTAL 

The films examined in this study were grown epitaxially by pulsed laser deposition on NdGaOa substrates, as 
detailed elsewhere [E3[. Each 1 cm x 1 cm film consists of buffer layers of semiconducting Pro.6Yo.4Ba2Cu30y_5 
that are 12 unit cells thick above, and 8 unit cells thick below, the YBCO film. The underlayer lessens the strain 
of substrate lattice mismatch on the YBCO, while the capping layer protects the ultrathin YBCO film from damage 
during handling. The three films we examined are nominally 4, 8, and 10 unit cells thick; however, it is likely that 
the top and bottom unit-cell layers are not perfectly smooth or homogeneous, leading to some uncertainty in the film 
thickness d. We therefore estimate d = 4 ± 1, 8 ± 1, and 10 ± 1 unit cells for the three films. 

We measure AJ X (T) with a two-coil mutual inductance technique described in detail elsewhere p4fl. The films are 
centered between two counterwound coils approximately 2 mm in diameter and 1 mm long. In a typical measurement, 



the sample is cooled to liquid helium temperatures, and a current of 100 [iA is driven at 50 kHz through the coil 
pressed against the back of the substrate. The supercurrents created in the sample are very nearly uniform through 
the film thickness. Thus, we observe the T c 's of all layers in the film. As the temperature slowly increases, the 
voltage induced in the secondary coil, which is pressed against the film, is measured continuously. The complex 
sheet conductivity cr(u>, T)d — ai(u), T)d — 102 (w, T)d of the entire film is deduced directly from the measured mutual 
inductance. The 2D penetration depth is obtained with an accuracy of about 3% from the imaginary part of the sheet 
conductivity from A J (T) = fiouia2(T)d, where /io is the permeability of vacuum. The uncertainty in film thickness 
enters only into calculation of the 3D penetration depth A _2 (T), and into the 2D penetration depth of each unit-cell 
layer. 

Figure 1 displays A^ 1 (T) as measured at 50 kHz for each film at optimal doping. An accelerated decrease in AJ 1 at 
high temperatures is most obvious for the 4 unit cell film, but similar features are apparent very close to T c for the 8 
and 10 unit-cell films as well. The transitions are slightly broader than in o-MoGe films B, presumably due to slight 
inhomogeneity. The vertical dotted lines mark the location of the peak in <7i(50 kHz,T), whose temperature we define 
to be T c and which very nearly coincides with the predicted KTB transition temperature, Tktb, for an isotropic 
film. To further motivate association of the drop in A J (T) with the KTB transition, we appeal to the frequency 
dependences of o\ and Aj 1 near the transition. Without going into detail, the observed shift of <j\ and Aj 1 to higher 
temperature with increasing frequency, shown in Fig. 2 for the 10 unit cell film, is quite similar to behavior seen in 
a-MoGe films [£2|. The experimental T c in these ultrathin films marks a KTB transition. 

Arrows in Fig. 1 indicate T KTB - the KTB transition temperature for uncoupled layers - showing that T KTB lies 
8 to 20 K below T c , a difference significantly larger than seen in YBCO crystals [7|n9| and simulations of J J arrays M. 
Another indicator of fluctuation effects is a rapid increase in downward curvature of A J for T > T KTB as predicted 
by the TPF model. The thin solid curves in Fig. 1 are extrapolations of quadratic fits to the data over T^ TB — 10 
K < T < T^tb- We emphasize that the quadratics fit the data extremely well over at least 20 K below T^ TB , 
highlighting the constant curvature of the data below T KTB . The insets to Fig. 1 show that only as T gets close 
to T c is a change in curvature of the data apparent. The proximity of the accelerated drops in A J to the predicted 
Tktb, and the lack of fluctuation effects (e.g., strong downward curvature) at T KTBl indicate the presence of strong 
interlayer coupling. 

We now address an important quantitative issue: the suppression of A J below its mean-field value just below 
Tktb- Taking the thin solid curves in Fig. 1 as reasonable estimates of the mean-field behavior of A^ 1 , we conclude 
that AJ 1 is 35% to 75% of its mean-field value at the transition, in agreement with results from a-MoGe films and 
calculations for 2D J J arrays M. 

Thus, the fluctuation behavior of ultrathin YBCO films is quantitatively similar to that of a-MoGe films. Within 
the TPF model, the only explanation is that interlayer coupling in the YBCO films is much stronger than expected 
from transport anisotropy. The strong a&-plane vs. c-axis anisotropy in transport coefficients could be explained 
phcnomenologically by occasional weak links between adjacent Cu02 bilayers, with most bilayers coupled strongly to 
their neighbors. 

III. CONCLUSION 

We have presented high precision, low frequency measurements of A^ X (T) and ct\(T) in YBCO films of 4 to 10 
unit cells thickness. The data indicate that interlayer coupling is strong, insofar as thermal phase fluctuations are 
concerned. Whatever the origin of the coupling, it is the presence of an additional layer or layers that lessens the 
effect of fluctuations on a single unit-cell layer. These results, and literature results for the aft-plane superfluid density 
of BSCCO crystals, argue that fluctuations are much weaker in cuprates than one would expect based on anisotropic 
transport coefficients. The power-law behavior near T c seen in the superfluid density of very clean YBCO crystals - 
cited as evidence of strong thermal phase fluctuations - is due to some other physics, perhaps competition between 
the superconducting order parameter and a normal-state order parameter that is sensitive to disorder. The idea of 
a normal-state order parameter has been proposed on various grounds Pq,M ■ The theory of Rokhsar H| suggests 
that the critical region may be more complicated than in the usual superconducting-to-normal transition when the 
normal and superconducting order parameters interact. 
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FIG. 1. XJ 1 vs. T measured at 50 kHz for three ultrathin YBCO films at optimal doping (thick solid lines). The vertical 
dotted lines locate the positions of the peaks in <ri(50 kHz,T), which very nearly coincide with Tktb- The thin solid lines are 
extrapolated quadratic fits to \~^ (T) over T K tb — 10 K < T < T^ TB and provide an estimate of the "mean-field" behavior of 
A^ 1 for T > Tktb- Insets are enlarged views of the transition regions. 
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FIG. 2. T dependence of the 2D penetration depth, \ ± , and conductivity, fiocodai, at 5 kHz (dotted lines) and 75 kHz (solid 
lines) for the 10 unit-cell-thick YBCO film. The arrow indicates T^ TB . 



